Quasiclassical and ultraquantum decay of superfluid turbulence 
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We address the question which, after a decade-long discussion, still remains open: what is the 
nature of the ultraquantum regime of decay of quantum turbulence? The model developed in this 
work reproduces both the ultraquantum and the quasiclassical decay regimes and explains their 
hydrodynamical natures. In the case where turbulence is generated by forcing at some intermediate 
lengthscale, e.g. by the beam of vortex rings in the experiment of Walmsley and Golov [Phys. Rev. 
Lett. 100, 245301 (2008)], we explained the mechanisms of generation of both ultraquantum and 
quasiclassical regimes. We also found that the anisotropy of the beam is important for generating 
the large scale motion associated with the quasiclassical regime. 
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The existence of a macroscopic complex order param- 
eter in superfluid helium ( 4 He and 3 He) constrains the 
vorticity to vortex lines, each line carrying one quan- 
tum of circulation n. This is in sharp contrast to or- 
dinary fluids, where vorticity is continuous. An impor- 
tant question is how quantum turbulence compares to 
classical turbulence JM- Experiments in helium have re- 
vealed two regimes [2|-|4[ of turbulent decay characterized 
by L ~ tr 1 (ultraquantum) and L ~ i~ 3 / 2 (quasiclassi- 
cal) behaviour, where t is time and the vortex line density 
(vortex length per unit volume) L measures the turbu- 
lence's intensity. In these two regimes the kinetic energy 
(per unit mass) decays as E ~ i -1 and E ~ tr 2 , re- 
spectively. (Here it seems appropriate to point to the 
detailed theoretical analysis |5| of energy decay in clas- 
sical, viscous, uniform and isotropic turbulence.) The 
second regime is thought to be associated with the clas- 
sical Kolmogorov distribution of kinetic energy over the 
length scales, but the nature of the first regime is still 
a mystery. Here we show that the first regime, associ- 
ated entirely with the Kelvin wave cascade along indi- 
vidual vortex lines, takes place when the energy input 
at some intermediate lengthscale is insufficient to induce 
the large-scale motion which is associated with quasiclas- 
sical, "Kolmogorov" turbulence. In other words, the first 
regime is a transient turbulent state which decays before 
energy can be transferred to large scales by vortex re- 
conncctions, which play a key role in this reverse energy 
transfer. 

Theoretical and experimental studies have revealed 
analogies between superfluid turbulence and classical tur- 
bulence, notably the same Kolmogorov energy spectrum 
in continually forced turbulence |6l-tl0(. as well as many 
dissimilarities and new effects. Our concern is the de- 
cay of pure superfluid turbulence at temperatures small 
enough that thermal excitations are negligible; in the ab- 
sence of viscous forces, in 4 He the only mechanism [l[ 
to dissipate kinetic energy is phonon emission at length 
scales much shorter than the average intcrvortex distance 



t L~ 1 / 2 . (In 3 He-B, which is a fermionic superfluid, the 
dissipation is thought to be associated with the Caroli- 
Matricon mechanism 11| of energy loss from short Kelvin 
waves into the quasiparticle bound states.) In this limit 
turbulence reduces to a very simple form: a disordered 
tangle of vortex lines, all of the same strength, moving 
in a fluid without viscosity, but still retains the crucial 
features of classical turbulence, the nonlinearities of the 
Euler equations and the huge number of length scales 
which are excited. 

By injecting negative ions in superfluid He in this 
zero-temperature limit, Walmsley and Golov ob- 
served two regimes of turbulence decay corresponding to 
two regimes of quantum turbulence discussed earlier in 



Refs. [12L Il3|, The negative ions (electron bubbles) gen- 



erated vortex rings 14J; the rings interacted with each 



other, forming a turbulent vortex tangle, which, in the 
first regime, decayed as L ~ t -1 . The second regime, 
characterized by L ~ i~ 3 / 2 , was observed if the injection 
time was longer. The same t~ 3 / 2 time dependence was 
observed in the spin-down of a vortex lattice fl5j ). and, 
at higher temperatures, during the decay of turbulence 
initially generated by a towed grid fl6j . Recently it has 
been also modeled numerically [17j ■ 

Walmsley and Golov argued that the second regime 
(which they referred to as Kolmogorov or quasiclassi- 



cal turbulence [12|, [13() is associated with the classi- 
cal Kolmogorov spectrum E% ~ fc~ 5 / 3 at wavenumbers 
k <C l/£, whereas the first regime (called Vinen or ul- 



traquantum turbulence [l2, 13 1) depends on energy con- 
tained at smaller scales, k ^> lji. The ultraquantum 
(L ~ t^ 1 ) and quasiclassical (L ~ t~ 3 / 2 ) decay regimes 
were also observed in 3 He-B by Bradley et al. [3j; in this 
case the turbulence was generated by a vibrating grid 
which sheds vortex loops in alternating directions. 

Following Schwarz (le| . we have modeled vortex lines 
as space curves s = s(t, £) (where £ is arclength) which 
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move according to the classical Biot-Savart law 



ds 

dt 
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where the line integral extends to the entire vortex config- 
uration C. Our model includes vortex rcconncctions and 



sound emission (for details see Refs |2l|-|23j]). The com- 
putational domain is a periodic box of size D — 0.03 cm. 
Modeling the experiments [H , the initial condition repre- 
sents the beam of vortex rings of radius R = 6 x 10~ 4 cm 
injected up to time t = 0.1 s with initial velocity ran- 
domly confined within a 7r/10 angle. 

The numerical techniques to de-singularize the Biot- 
Savart integral, discrctize the vortex filaments over a 
large variable number of vortex points Sj, and per- 
form vortex reconnections are standard in the litera- 
ture fl8j | . Details of our algorithms are in our previ- 
ous papers [21 - 23 , which also describe the tree algo- 
rithm used to speed up the calculation of Biot-Savart in- 
tegrals. Our model includes small energy losses at vortex 
reconnections, as described by more microscopic calcula- 
tions based on the Gross-Pitaevskii equation [24j . Energy 
losses due to sound emission are modelled by the spatial 
discretization [2l| : vortex points are removed if the local 
wavelength is smaller than a given minimum resolution 
5 = 5 x 10~ 4 cm. In our calculations, the time step is 
At = 5 x 10~ 6 s. We have tested that the ultraquan- 
tum and quasiclassical behaviours remain the same if S 
is halved. 

By numerically integrating Eq. (|T|) we have found that 
the vortex rings interact, reconnect and, as envisaged by 
Bradley et al. (20| . form a tangle; the vortex line density 
L reaches a peak and then decays, see Fig. Q] (left), in 
agreement with the observed ultraquantum (L ~ t -1 ) 
behaviour. During the decay, the kinetic energy (per 
unit mass), E, has the expected E ~ t^ 1 behaviour. 

We have repeated the calculation with longer injection 
time, up to t = 1 s. The peak value of L is thus about 
10 times larger than in the ultraquantum case, as in the 
experiment [2j. We have found that, as shown in Fig. Q] 
(right), after the initial transient the decay assumes the 
quasiclassical (L ~ r, -3 / 2 ) form observed in the exper- 
iments 0, [3J]. We have also checked that E ~ t~ 2 , as 
expected. The same quasiclassical and ultraquantum de- 
cays are obtained with half the numerical resolution along 
the vortex filaments. 

It should be emphasized that left and right of Fig. [1] 
do not represent different stages of turbulence but repro- 
duce two different experiments resulting, respectively, 
in two different regimes of decay: ultraquantum and qua- 
siclassical. The key parameter, determining which of the 
two regimes will be realized, is the time of injection of 
vortex rings. 

Assuming the classical expression dE/dt = —vuj 2 , 
where w is the vorticity, and the identification lo = kL, 
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FIG. 1: Vortex line density L (cm -2 ) vs time t (s) cor- 
responding to small (left) and large (right) times of vortex 
rings injection. Fitting the decay as L ~ t~ a , we obtain 
a — 1.07 and 1.48, respetively. The dashed lines show the 
ultraquantum L ~ t" 1 (left) and quasiclassical L ~ t -3 / 2 
(right) behaviours. 



we interpret the results in terms of an effective kine- 
matic viscosity v, which we call vy ("Vinen") and vk 
( "Kolmogorov" ) respectively for the two regimes 0]. 
The values of the effective kinematic viscosities vy and 
i>k have been obtained as in Ref. Q by fitting respec- 
tively L w B/(vyt), where B = (1/(4tt)) In (£/a ) and 
cio « 10~ 8 cm is the vortex core radius, and L w 
($C) 2/ ' i v~ K 1/2 K- 1 ki 1 t-' i / 2 , where 27r/fci is the large scale 
and C = 1.5 is the Kolmogorov constant. In applying 
these formulae we have taken into account the fact that 
for the calculations presented here in the ultraquantum 
case the computational box is not entirely full, and that 
in the quasiclassical case the largest length scale is of 
the order of 0.06 cm, as visible in PDF(C). We obtain 
vy I k pa 10 -1 and vk/k « 10~ 3 , which compare fairly 
well with Walmsley & Golov's vy/n « 0.08 to 0.1 and 
v k /k w 0.002 to 0.01. 

To understand the nature of the two regimes we have 
examined the time behaviour of the probability density 
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function PDF(C) (normalized histogram) of the local 
vortex line curvature, C = \d 2 s/d£ 2 \. In both ultraquan- 
tum and quasiclassical case, the initial PDF develops in 
time to larger and smaller values of C . In the quasiclas- 
sical case, however, there is a much greater build up at 
small values of C, see Fig. [3J this means that, as the ini- 
tial vortex rings entangle, large-scale structures are cre- 
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FIG. 2: Probability distribution functions of the vortex line 
curvature, C (cm - ). Solid line: ultraquantum decay at t — 
0.08 s; dashed line: quasiclassical decay at t — 0.7 s. 

ated consisting of long vortex filaments which can extend 
across the entire computational domain. 

This generation of large length scales is apparent in 
Fig. [3l where we show the evolution of the kinetic energy 
spectrum Ek, defined by 

1 f 1 f°° 

E= V Jy2^ 2dV = J Ekdk ^ 

(where V is volume and k the magnitude of the three- 
dimensional wavevector k). In both ultraquantum and 
quasiclassical cases the energy is initially concentrated 
at intermediate wavenumbers. It is apparent (see Fig. [3] 
right) that in the ultraquantum case the value of k = k* , 
where Ek has the maximum does not change, and the 
ratio of energy transferred to large scales, J^* Ek dk, to 
that transferred to small scales, J k °° Ek dk, remains small 
(< 0.13) at all times. In the quasiclassical case, how- 
ever, a significant amount of energy is transferred to 
small wavenumbers, leading to the formation of the Kol- 
mogorov fc~ 5 / 3 spectrum, see Fig. [3] The spectrum main- 
tains the Kolmogorov scaling during the decay stage, con- 
sistently with the observation that L ~ i -3 / 2 , even for 
relatively small values of L which would otherwise decay 
as i -1 if L were small initially. 

To interpret these results we remark that in both ex- 
periments 0, 0] the initial vortex rings do not move 
isotropically, but essentially travel in the same direction 
as a beam. This anisotropy is important in creating 
large length scales, provided that the initial density of 
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FIG. 3: (Color online) Energy spectrum (arbitrary units) vs 
wavenumber k, cm -1 corresponding to the quasiclassical de- 
cay (left) with the solid line at t — 0.1 s, dot-dashed t = 1.1 s, 
and dashed t = 3.4 s, and to the ultraquantum decay (right) 
with the solid line at t — 0.07 s, dashed t = 0.12 s, and 
dot-dashed t = 0.6 s. In the former, note the formation and 
decay of the Kolmogorov spectrum (indicated by the straight 
dashed line). 



the rings is large enough. The argument is the follow- 
ing. Energy and speed of a vortex ring of radius R are 
respectively proportional and inversely proportional to 
R. Consider the collision of two vortex rings of approx- 
imately the same size. If the collision is head-on, the 
outcome of the reconnection will be two vortex loops of 
approximately the same size, as shown schematically in 
Fig. [4] (bottom). If the two rings travel approximately 
in the same direction, the reconnection will create two 
vortex loops, one small and one big, see Fig. H] (top). To 
test the idea that an anisotropic beam facilitates the cre- 
ation of length scales, we have performed numerical cal- 
culations in which the initial distribution of vortex rings 
differs only by the orientation of the rings: in one case 
the rings pointed isotropically in all directions, and in an- 
other case they pointed in the same direction. Figured] 
confirms that the anisotropic initial condition generates 
smaller values of curvature (that is, larger length scales) 
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FIG. 4: Head-on reconnection of two vortex rings of sim- 
ilar size tends to produce two vortex loops of similar size 
(bottom), whereas the reconnection of two rings traveling in 
the same direction produces one vortex loop which is much 
smaller and one which is much bigger (top). 




10 4 10 3 10 2 10 ' 10° 

R 



FIG. 5: Probability distribution functions, PDF(Ti) of loop 
sizes R, (cm). Dashed line: initial PDF. Grey line: resulting 
PDF for isotropic initial condition: no large loops are cre- 
ated by the vortex reconnections. Black line: resulting PDF 
for anisotropic initial condition: note the generation of large 
loops. 



as well as bigger ones. To highlight the geometrical role 
of vortex reconnections, the calculation was performed 
by replacing Eq. (JTJ with the local induction approxi- 
mation ds/dt = /3s' x s" (where the prime denotes 
derivative with respect to arclcngth and /3 is constant); 
in this way the rings interacted only when they collided. 
The result was similar to that obtained using the full 
Biot-Savart calculation. 

In conclusion, our model reproduces both the ultra- 
quantum [L ~ t~ x ) and the quasiclassical (L ~ £~ 3 / 2 ) 
turbulent decay regimes which have been observed in 
the experiments [H, Q and explains their hydrodynam- 
ical natures. By examining the curvature PDF and the 
energy spectrum, we have found that in the quasiclassical 
regime the initial energy distribution is shifted to large 
scales, and a Kolmogorov spectrum is formed. In the ul- 
traquantum case, the spectrum decays without this en- 
ergy transfer. In the case where turbulence is generated 
by forcing in the vicinity of some (intermediate) length- 



scale, as in the experiments 0, H|, we found that the 
ultraquantum regime is induced only if the total energy 
input is relatively low, while the higher energy input (by 
e.g. theprolonged injection of the vortex rings in exper- 
iments 0,3) generates the large scale motion and hence 
the quasiclassical, Kolmogorov regime of turbulence. We 
have also found that the anisotropy of the beam of vortex 
rings is important, as reconnections of vortex loops trav- 
eling in the same direction are very effective in creating 
larger length scales. 
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